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Abstract 

A dynamical model based on effective Lagrangians is proposed to describe the bremsstrahlung 
reaction ttN — > irN'-f at low energies. The A(1232) degrees of freedom are incorporated in a way 
consistent with both, electromagnetic gauge invariance and invariance under contact transformations. 
The model also includes the initial and final state rescattering of hadrons via a T-matrix with off 
the momentum-shell effects. The double differential distribution of photons is computed for three 
different T-matrix models and the results are compared with the soft photon approximation, and 
with experimental data. The aim of this analysis is to test the off-shell behaviour of the different 
T-matrices under consideration. Finally an alternative simpler dynamical model that incorporates 
the unstable character of the isobar-A(1232) through a complex mass, is presented. As we will see 
it is suitable for the study of the magnetic moment of the resonance. 
PACS numbers: 25.80.Ek, 13.60.-n, 13.75.-r 
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1. INTRODUCTION 



In order to extract resonant parameters of the nucleon resonances(N*) from the ^yN — > ttN reaction, 
it is important to evaluate the background contribution to isolate the resonant peak. An important 
contribution to the background of this photoproduction reaction is provided by the final state rescat- 
tering (FSI) of the nN system |], @, ||, |J • Consequently, we require the knowledge of the T- matrix 
that describes this rescattering process, in the off the momentum-shell regime(off-shell),i.e, we need 
T(q <f; z(q)) with \q\ ^ \q '\ where z is the total energy of the nN system as a function of the 
relative momentum |g| of the initial state. This particular rescattering amplitude is more properly 
called half-off-shell T-matrix. In all cases the so called 'realistic' interactions are fitted to reproduce 
the phase shifts in elastic nN scattering, which only depends on the on-shell (\q\ = \q f \) values of 
the relative momenta. Thus, elastic scattering is not useful to constrain the off-shell behavior of the 
T-matrix because interactions working similarly in elastic scattering, may have different behaviour 
in the off-shell regime. 

Another reaction where the irN off-shell T-matrix is required is ttN — > nN'-f bremsstrahlung. 
This process has been studied within the Soft Photon Approximation (SPA)||. Within this approx- 
imation the full amplitude, expressed as a power expansion in the photon energy, depends only on 
the electromagnetic static multipoles of tt and N and on the T-matrix (and its derivatives) of the 
corresponding non-radiative process ||. The SPA reproduces very well the experimental data on 
radiative ttN scattering, in spite of some objections that have been raised recently |7| regarding the 
departures of the formulation of ref. H from the original Low's prescription. Because of the power 
expansion in the photon energy, the total SPA amplitude depends on the derivatives of the T-matrix, 
and thus on off-shell effects. Since the non-leading terms (of the power expansion) are fixed by im- 
posing electromagnetic gauge invariance, the off-shell effects of the T-matrix cancel up mutually. 
Therefore, the information on the off-shell behaviour of the T-matrix can be tested by adding the 
contributions to the radiative nN scattering within the framework of an specific dynamical model. 

The purpose of the present paper is to check the off-shell behaviour of three different T-matrices for 
ttN rescattering in the reaction ttN — ► 7riV7. We use a dynamical model to describe the ttN — > nN'-f 
reaction. The gauge invariant electromagnetic current is constructed explicitly, with vertices and 
propagators derived from the relevant hadronic and electromagnetic Lagrangians. We include also 
two-body meson exchange currents, and the full energy-momentum dependence of the T-matrix for 
the elastic ttN scattering which exhibits its off-shell behaviour. Finally we implement this model 
with different T-matrices in order to compare their different off-shell dependence. 

On the other hand since different T-matrices depend on many parameters (bare masses and 
coupling constants of the hadrons, cut-off form factor parameters, etc.), results difficult to use a 
dynamical model based on a T-matrix input for analyzing resonance unknown properties. For this 
reason we also present a simpler formalism that assumes the isobar-A as the main degree of freedom, 
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and gives it an unstable character through a complex mass. This approach will be used to study the 
anomalous magnetic moment of the resonance. 

This paper is organized as follows. In section 2 we will construct the gauge-invariant amplitude 
of radiative irN scattering. The Lagrangians used to construct the gauge invariant current of our 
process, are provided also in this section. The second simpler dynamical model will be described in 
section 3 . Finally the results and conclusions are given in section 4. 



2. GAUGE INVARIANT BREMSSTRAHLUNG AMPLITUDE AND 

DYNAMICAL MODEL 



In the pion-nucleon bremsstrahlung process we deal with a problem of scattering by two potentials 
||: the strong pion-nucleon and the electromagnetic interactions. The cross section for ttN — > ttN'j 
process reads 



da 



dk 



dqf f dpf 



^ 7 J J E (Pf) 



(27r) 4 5 4 (p i + q i -p f -q f -k) 



ex,ms f ,msi 



m 



N 



M 7rN7t7rN (ex, k;qf,pf,msf,qi,pi, msi) 



where q = (00, q), p = (E,p) and k = {u^^k) denote pion, nucleon and photon four-momenta, 
respectively; ms is the nucleon's spin projection and e\ indicates the polarization four-vector of the 
photon. The subindexes i, f refer to initial and final state quantities. 
The Lorentz invariant amplitude^ M n N 1<n N explicitly reads 



M nNl)nN = {u{pf,msf)\M^ N ^ N {e x ,k]qf,pf,q h Pi)\u{p h mSi)), 



(2) 



where u(p,ms) denote nucleon Dirac spinors, and the amplitude operator M^n^n is obtained from 
the coupled channel Bethe-Salpeter equation for the tiNj system as follows (we consider electromag- 
netic interactions at the lowest order) 



M, 



ttNj,ttN 



+i 
+i 



V-ttN'y^N 

dq A 



2 j dq 4 dq 14 



(2vr) 4 (2vr 



V-KN-r^N^G^N^M^N^jy^q) + M 1TN ^ N {q)G 7T N{q)y-KN-y,TTN{q) 

M n N,nN (q ) Gtj-jv (q ) VirNy,irN (q , q) G-kN (q)M 7r N,TrN(q) 



(3) 



throughout this paper, M will denote the amplitude generated by the operator M,ie., M — (u\M\ 



u . 
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In terms of the above operator amplitude the T-matrix, defined as 



T(qf,Pf,qi,Pi) = 7^pM7rJV,7riv(g/,P/;gi,Pi), (4) 



satisfies the integral equation 



U + i/ ^U{q)G{q)T{ql (5) 



G = (2ti) 3 G 7tN . 

In the previous equations Vy denote M-matrix elements corresponding to the irreducible Feynman 
diagrams for each process, while Gi is the product of Feynman propagators of intermediate particles. 

Following the Thompson's prescription |J, we can set the above integrals in a three-dimensional 
form as follows (we set in the center of mass frame of the ttN system) 

f(q',q,z) = U(q',q) + J d 3 q n U(q',q")GTH(z,r)f(q\q } z), (6) 

with, 

mN E \u(-q n ,ms n ))(u(-q n ,rns n )\ 
° ThM) = 2u{q^)E{-r)~ z-z»+i V ' (?) 

where z" = E(—q") + uj(q "). In the above expressions Gth denotes the Thompson propagator 
replacing the full G n N Feynman propagator which, as a consequence of the three-dimensional reduc- 
tion, eliminates the propagation of antiparticles and put intermediate particles on their mass-shell. 
The kernel function U(q q) contains all the 7r./V-interaction irreducible diagrams to be iterated in 
the T-matrix calculation, but usually only second-order contributions are kept. 
The electromagnetic current V^n^^n can be broken into two pieces 

Kr7V7,7riV = KrJV 7)7r JV + ^kN^itN) (^) 

where the upper indices denote one- and two-body contributions, respectively which are obtained by 
coupling the photon to all the internal lines of U. As is known, the operator t^v 7)7rJV mus t be added 
to the electromagnetic current in order to satisfy the electromagnetic gauge invariance of the total 
amplitude |T0[|, while the one-body amplitude Vj^ y ^ wN vanishes for free hadrons. Both contributions 



to the total amplitude are illustrated in Fig. 1. 
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Fig.l One- and two-body contributions to the bremsstrahlung current amplitude. 



We follow a procedure that put the bremsstrahlung amplitude manifestly gauge invariant. By 
replacing eqs. and (||) into eq.(|3]), only in the one-body component of the amplitude M^f %7rN = 

MnN-y^NiV^^^x), we can isolate the lowest order nonzero contribution of the one-body current. 
After the three dimensional reduction, the total amplitude can be rewritten as 



with 

Kn-y,*n = (u(-q ' - k/2, ms f ) \V^nj,tvn{^x, k, q ', q) \u(-q, msj) , 

Kn\,n = J dq m (u(-q\ms f )\f^(q\?,z')G TH (z\^^ 

M^^n = J dq" 3 (u(-q' - ms/) (^^^(eA, fc, g ^ 9, ^) mSi)), 

M£%Z N = J dq" 3 J dq"' 3 (u(-q' - k/2,ms f )\f^(-q' - k/2,$',z') 

G TH (z', q")VnN 7 ,-KN(e\, k, q", q")G T H(z, q")f(q", q, z)\u(-q, msi)) (10) 
where the current 

VVa^ttTv = iV}lf linN GU + i&GV^^N + V^^ N , (11) 

generates a gauge invariant Born amplitude V^N-y^N, that involve all the possible ways of attaching a 
photon to the ttN scattering amplitude U, and contains the full propagator G ~ G^n- The operator 
T^(z'), where z' = z + cu 7 , obeys eq.@ if we change r\ — > —i] in eq. (|7|). 

The superscript pre (post) in Eq. (9) indicate that the photon is emitted before (after) the action 
of the T-matrix, while the superscript double refers to a double-scattering term where the photon 
is emitted from internal lines between two T-matrices. In the above equations the Born, pre and 
double amplitudes were evaluated in the initial center of mass frame (g) = q ' — k/2,pf = —q' — k/2, 
iji = —pi = q), while the other (post) amplitude was evaluated in the corresponding final frame((f/ = 
—pf = q ', (ji = q + k/2, Pi = —q + k/2). The different terms in Eq. (9) are illustrated in Fig. 2a 
(KrV7,7rAf) and in Fig. 2b (remaining terms). 

Fig. 2 (a) Gauge- invariant bremsstrahlung current amplitude, (b) Post-, pre- and double-scattering 
amplitude contributions (se eq. (9)). 
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Within the SPA, the total bremsstrahlung amplitude can be split into external (E) and internal 
(J) contributions: 

M nNl ^ N = M^ 7i7rJV + M wNy ^ N , (12) 

where we can identify 

and the internal contribution M^ Nry ^ nN can be obtained "by enforcing" the gauge invariance condition 

Mkn^nW = fe") = 0. (13) 

The SPA the M^n^n amplitude depends only on the elastic T-matrix, because derivative terms 
of T cancels in the addition of internal and external contributions. Let us emphasize that any 
dependence on the structure of internal contributions (in particular, the dependence of off-shell 
effects) are of higher order in cj 7 and must be included explicitly in the amplitude in a gauge invariant 
way 

The bremsstrahlung amplitude will be computed along the lines developed in eqs. (9-11), using 
a potential U obtained from effective Lagrangians, and three specific models for the T-matrix to 
describe the irN rescattering. 

The three models advocated for the T-matrix will be called OBQA, SEP and NEW, respectively. 
The OBQA version for the T-matrix interaction [jllj, is based on a model that includes it and p 



mesons exchange through a correlated 2n exchange potential. The SEP model for ttN rescattering is 
generated [12|] by a phenomenological separable potential. Finally, the NEW model is obtained 



from exchange of it and (sharp) p mesons. The operator U is constructed from a Lagrangian density 
that includes the nucleon(iV), the A-isobar, and the it, p and a mesons 



crNN (T7T7T J 

(14) 

while the electromagnetic currents are obtained from the hadronic Lagrangian density through min- 
imal coupling of the photon. The Scattering amplitude U is depicted in Fig. 3, while the amplitude 
Kr7v 7 ,7r7v can be obtained by coupling the photon to all diagrams in U as shown in fig. 4. 

Fig. 3 Born amplitude corresponding to the 7rN potential operator U. The first diagram denote 
the nucleon-pole, and the A-pole corresponds to the third diagram. The fifth and sixth diagrams 
correspond to p and a mesons exchange. 

Fig. 4 The gauge-invariant amplitude obtained by coupling a photon to the Born terms in fig. 3, 
together with the two-meson exchange currents (fifth, sixth and seventh graphs in each line), (a) 
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Contributions obtained from the N and A direct-pole diagrams in fig. 3; (b) Terms generated by the 
cross-pole diagrams in fig. 3, and (c) diagrams obtained from p and a exchange contributions in fig.3. 



The good convergence properties of the scattering equations given in Eqs. fllOD can be obtained by 
introducing hadronic form factors, which are supposed to describe the composite nature of hadrons. 
It is a common practice to use different parametrizations of the form factors for different T-matrices. 
For example, the OBQA model uses monopole and dipole forms with cutoff parameters ranging from 
A = 1200 — 1600Mey ||ll|| . In the case of the SEP interaction different form factors are introduced 
for each partial wave component |12|], while in the NEW form factors usually advocated are of 



monopolar form with A = 1300 — 2300M eV However, the introduction of form factors replacing 
point vertices in V^n^^n spoils the gauge invariance of the total amplitude. Fortunately, the gauge 
invariance of the amplitude can be recovered by using the method of Gross and Riskafl5| which, 
however, does not yield unique electromagnetic couplings to hadrons. Therefore, we follow the more 
simple prescription of using a common form factor jl6[ of monopole type 



A 2 



A 2 + <f 2 ' 

where the scale A can be adjusted at a given incident energy for each T-matrix model. 



(15) 



3. DYNAMICAL MODEL II 



The dynamical model described in the previous section has the adventage of decribing irN -FSI 
with great detail by the inclussion of many components in the T-matrix potential. Nevertheless, 
these various effective contributions depend on time on certain parameters that are fitted in order to 
reproduce scattering phase shifts. This makes difficult to analyze (by fitting) unknown parameters of 
the resonances using a scattering matrix, on which we have additional parameters to be controlled. 
In this section we will adopt an alternative dynamical model that is simpler, and gives us a tool to 
study the, until this moment, unknown magnetic moment of the isobar-A(1232) resonance. 

Taking into account that in the range of energies we will consider (7] a 6 ~ 300MeV) the main 
contribution comes from the A-intermediate terms, since we are in the resonance region, we only 
consider the A-pole contribution in the hadronic potential (third contribution in fig. 3). This reads 

(j = f (vriV -> A)^f(A -> ttA^, (16) 

where 

G^(q) = 1 2 6^{q) (17) 
q z — m% 

7 



0^(q) = U+m A ) 
2 



3m 2 A 



1 2 q^q v 1 q^Y - q v Y 

6 3 m A 6 m A 

(q 2 - m 2 A ) \ (7V - tV) + {4 + m A )Yl 



is the A propagator jL% . 

Iteration of this contribution to all orders through the T- matrix eq. (0) leads to 



is) 



U = f (vriV -> A) v g u ^t(A -> vriV), 



where 



(19) 



r d°*"( ff ) 



- m A )<?£ - E( ff )g' 
is a modified A propagator with a self-energy 

S(ff)S = &*{q) J d 3 pT(A -> 7riV)(p)^(p,Z)f(7riV 



A)(p)c 



(20) 



In order to simplify we adopt the renormalization recipe 

mlgZ + Z(qr a K(M 2 A -iM A r)gZ 

with Ma and T constants. This should be equivalent to make the replacement m\ — > M\ — iM A T 
in the denominator of the G^lq), with which we are considering the unstable character of the A. In 
order to get gauge invariance at the moment of coupling the photon to each line of the potential, we 
must make this replacement in the full propagator. Finally to calculate the bremsstrahlung amplitude 
within this model we must replace U by U in eq. QTTl) and keep only the born contribution in eq.(§). 

The anomalous magnetic moment of the A, which is not wellknown will be moved in orther to 
fit the region of high-energy photons since in the A electromagnetic vertex 



VfiOL 



1 1 1 



2m, 



(21) 



where k A and e A are the anomalous magnetic moment^ and charge operators whose action upon the 
Rarita-Schwinger field give as eigenvalues the corresponding values of these properties, the k A term 
goes as c<j 7 . 



2 We restrict ourselves to the A ++ contribution, the only one for which we have experimental information on ka 
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4. NUMERICAL RESULTS AND CONCLUSIONS 



The differential cross section da / ' dVt 1T dVt 1 duJ 1 to be compared with experimental data can be 
obtained from eq.(|I|), where the amplitude M^N-y^N is calculated from eqs.(9-ll). The dynamical 
model approximation (DMA I) advocated in the present paper contains the following steps. The 

current operator V n ^ 7:7T N is computed from the effective Lagrangian described in section 2, and the 
monopole form factor given in eq.(^) to have good convergence of the intermediate momentum 
integrals. 

As is known, double scattering terms have a significant contribution in the case of proton— proton 



bremsstrahlung, mainly in the end point region of co> 7 ||19fl . In the present work we will neglect M^^ N 



in eq. (|T0|) ) , because the numerical calculation of the two three-dimensional integrals requires an 
enormous computation effort. Nevertheless, we keep double scattering-like contributions in post and 
pre amplitudes coming from current components iU^GV^ y7rN , and iV^j^^GU, respectively. In 
order to compare the approaches provided by the DMA I and SPA, we fix M^n^^n to coincide 
quantitatively at low photon energies. 

In addition we are going to study the sensibility of the bremsstrahlung cross section with the 
value of the A magnetic moment. For this purpose we also evaluate the cross section through the 
alternative dynamical model described in section 3 (DMA II). 

For illustration purposes, we will implemented the DMA I approach with the OBQA, SEP and 
NEW T-matrices for the specific example of n + p — > 7i + p-f. The different coupling constants and 
masses needed to evaluate V nN7i7rN were taken from the model II of ref. from fll2"| , and from 



I3| . For direct pole diagrams we use bare masses and coupling constants since they get dynamically 



dressed in the T-matrix scattering eq.(|5|). In the OBQA case we replace the 2ir correlated exchange 
potential, by 7r and p sharp mass exchange terms, since they do not lead to sizable differences as 
shown in ref. |nj. In the SEP case we use the same coupling constants and masses, because the 
scattering potential is not generated from a dynamical model. 

We will compare the theoretical predictions to the experimental cross sections measured by 
Nefkens [ [20fl ( EXP ) , which have been reported for different kinematical configurations. In EXP the 
pions were detected at fixed angles 9 n = 50.5°, <p n = 180°, for three different energies of incident 
pions (269, 298 and 324 MeV), and the photons were detected at various 6* 7 ,0 7 angles in the range 
of energies ^ = — 150 MeV. 

Our results for the cross sections in the DMA I approach, for the different ansatz of the T- 
matrix, are shown in Fig. 5. The SPA and experimental values are also plotted for comparison. 
The predictions of the DMA I approach shown in fig. 5 using the three models of the T-matrix, 
are compared to the results of EXP for photon angles given by Gu = # 7 = 103°, 7 = 180°. Since 
the parameters entering T-matrices are usually quoted to reproduce the elastic phase shifts, the 
cutoff parameters A were fixed in order to have good coincidence between the predictions of the 
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different interactions and the SPA at low c<j 7 values. We get (in MeV units) Aobqa = 750, 700, 600, 
A SEP = 700,600,500, and A NEW = 550,500,450, for incident pion energies T lab = 269,298,324 
respectively. Observe that the value of A found in the case of the SEP interaction for Ti ab = 298MeV 
is consistent with the one previously found for pion photo-production at Ti ab = 300MeV[]12|, while 



the OBQA values are roughly consistent with the form factors used in ref. JTTj] , which corresponds to 
a monopole form-factor with A w 800MeV[lf|. 



Fig. 5 7rA7 cross section for T lab = 269, 298 and 32AM eV and G u = ft, 



103° 



180° 



in 



EXP, calculated in the DMA I for the different T-matrices. 
and the measured values. 



We also include the SPA cross section 



Results for the cross section within the DMA II are shown in Fig. 6 for different photon detectors, 
G14 = # 7 = 103°, 7 = 180°;G n = 6^ = 160°, 7 = 180°; and G 7 = # 7 = 120°, 7 = 0°. We will fix 
M A = \232MeV and T = HOMeV as reported in previous experimental works fl2T||. The coupling 
constant /na-k, was determined fitting the low-energy photon emission region (elastic scattering). 
We obtained Jna-k = 0.24, which is in the range of previous works. Independently, the A magnetic 
moment ka was shifted in order to show how the high energy photon region can be fitted. 

Fig.6 7tA 7 cross section for T lab = 269MeV and G 7 = 6 1 = 120°, 7 = 0°, G n = Q^ = 160°, 7 = 
0°, and G?i4 = # 7 = 103°, 7 = 180° in EXP, calculated in the DMA II for two different values of k a . 

In almost all the cases, the predictions of the DMA I (fig. 5) lies above the experimental cross 
section and the SPA for energies ou 1 > 20MeV. One of the reasons for this may be the use of an 
overall form factor to cure the gauge invariance problems. The total bremsstrahlung amplitude is 
built up, as can be seen from eqs. @ and fllOD , by adding different components. It is not expected 
that the common form factor works satisfactorily by adding up these components as it does the one 
used to generate the individual T-matrices, which change their values from vertex to vertex. The 
comparison between the results of the SPA and DMA I schemes shows that the additional off-shell 
effects, added coherently to the lowest order contributions, may have important contributions since 
they do not cancel exactly with the derivate terms of the T-matrix appearing from the soft-photon 
expansion. 

From fig. 5 we can check that the SEP interaction provides the closest results to the experimental 
cross section with a departure starting for cu 7 > AOMeV. This indicates that the dynamical model 
involved in the SEP interaction gives the smallest off-shell effects. On the other hand, the strongest 
off-shell effects appear in the OBQA model. This conclusion agrees with a previous study 0] on 
observables in pion photo-production experiments. 

As it was discussed in section 2, the off-shell contributions to the external and internal amplitudes 
within the SPA cancel each other, thus we cannot study these effects within this approximation. In 
addition since we get gauge invariance in the SPA by adjusting the internal amplitude, the gauge- 
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invariant electromagnetic currents remain hidden. In the DMA I approach, these cancellation must 
occur explicitly between the different components of the amplitude (the so-called born, pre, post con- 
tributions). The departure of the different T- matrices from the SPA can be used to estimate the 
size of unbalanced off-shell terms and provide a test of their off-shell behaviors. Also, since the 
electromagnetic gauge-invariant current is constructed explicitly from effective Lagrangians, we can 
use the radiative irN reaction to study the relevance of the degrees of freedom and the parameters 
involved in this dynamical model. 

On the other hand from fig. 6 we can see as this reaction within the DMA II (tree-label + complex 
Ma), is sensible to the changes of ka- This model seems appropiated to analyze the anomalous 
magnetic moment of the A resonance, since within the SPA the the 'enforcing' gauge invariant 
procedure does not fix this term, which is self-gauge invariant. In order to give a realistic value , a 
more carefull fitting procedure of Ma, T, and k&, must be done. 
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